library(stargazer)
library(tidyverse)
library(xtable)
library(clubSandwich)
library(lmtest)
library(rddtools)
library(sandwich)
library(lmtest)
library(usmap)
library(rgdal)
library(tigris)
library(stringr)
library(tidycensus)
library(acs)
library(sf)
library(tmap)
library(interflex)

#load data

ag <- read.csv("data/ag.csv")
national_data <- read.csv("data/national_data.csv")
national_data_sheriff <- read.csv("data/national_data_sheriff.csv")
new_jersey_data <- read.csv("data/new_jersey_data.csv")

new_jersey_data$condition <- as.factor(new_jersey_data$condition)

new_jersey_data$condition <- factor(new_jersey_data$condition,levels(new_jersey_data$condition)[c(2,1,3,4)])
national_data$condition <- as.factor(national_data$condition)
#####################################
#####################################
### Main Manuscript
#####################################
#####################################

#####################################
### Figure 1 - Experiment Coverage
#####################################

source("figure1.r")

#####################################
### Figure 2 - Response Correlates
#####################################

source("figure2.r")

#####################################
### Figure 3 - Response by Treatment
#####################################

source("figure3.r")

#####################################
#####################################
### Appendix
#####################################
#####################################

#####################################
### Table C1 - Census and Agency Data
#####################################

# Generation of this table requires agency identifiers that cannot be anonymized. 
# Our IRB protocols mandate as much anonymity as possible for sampled agencies.  
# We have omitted this code to comply with our IRB ruling.  We note that this is
# a descriptive table that plays no role in our analysis or results.

#####################################
### Table D1 - NJ Balance
#####################################

source("tabled1.r")

#####################################
### Table D2 - National Balance
#####################################

source("tabled2.r")

#####################################
### Table E1 - Sampled Agencies by State
#####################################

source("tablee1.r")

#####################################
### Table E2 - Estimates from Figure 2
#####################################

source("tablee2.r")

#####################################
### Figure E1 - Alternate Figure 2
#####################################

source("figuree1.r")

#####################################
### Table E3 - Estimates from Figure 3 (NJ)
#####################################

source("tablee3.r")

#####################################
### Table E4 - Estimates from Figure 3 (National)
#####################################

source("tablee4.r")

#####################################
### Table E5 - Trump Vote Share and Response
#####################################

source("tablee5.r")

#####################################
### Figure E1 - Trump Vote Share and Response
#####################################

source("figuree2.r")

#####################################
### Figure G1 - Performance and Treatment
#####################################

source("figureg1.r")

#####################################
### Figure H1 - Sensitivity 
#####################################

source("figureh1.r")

#####################################
### Table J1 - Alternative Coding
#####################################

source("tablej1.r")

#####################################
### Table k1 - Chiefs v. Sheriff
#####################################

source("tablek1.r")
